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We present a detailed quantitative description of the recently proposed “slingshot effect”. Namely, 
we determine a broad range of conditions under which the impact of a very short and intense laser 
pulse normally onto a low-density plasma (or matter locally completely ionized into a plasma by the 
pulse) causes the expulsion of a bunch of surface electrons in the direction opposite to the one of 
propagation of the pulse, and the detailed, ready-for-experiments features of the expelled electrons 
(energy spectrum, collimation, etc). The effect is due to the combined actions of the ponderomotive 
force and the huge longitudinal field arising from charge separation. Our predictions are based on 
estimating 3D corrections to a simple, yet powerful plane 2-fluid magnetohydrodynamic (MHD) 
model where the equations to be solved are reduced to a system of Hamilton equations in one 
dimension (or a collection of) which become autonomous after the pulse has overcome the electrons. 
Experimental tests seem to be at hand. If confirmed by the latter, the effect would provide a new 
extraction and acceleration mechanism for electrons, alternative to traditional radio-frequency-based 
or Laser-Wake-Field ones. 


I. INTRODUCTION AND SET-UP 

Laser-driven Plasma-based Acceleration (LPA) mech¬ 
anisms were first conceived by Tajima and Dawson in 
1979 [1] and have been intensively studied since then. In 
particular, after the rapid development of chirped 

pulse amplification laser technology - making available 
compact sources of intense, high-power, ultrashort laser 
pulses - the Laser Wake Field Acceleration (LWFA) 
mechanism [D m E] allows to generate extremely high 
acceleration gradients (>lGV/cm) by plasma waves in¬ 
volving huge charge density variations. Since 2004 ex¬ 
periments have shown that LWFA in the socalled bub¬ 
ble (or blowout) regime can produce electron bunches of 
high quality (i.e. very good collimation and small en¬ 
ergy spread), energies of up to hundreds of MeVs [6]-[8] 
or more recently even GeVs laiTo]. This allows a rev¬ 
olution in acceleration techniques of charged particles, 
with a host of potential applications in research (parti¬ 
cle physics, materials science, structural biology, etc.) as 
well as applications in medicine, optycs, etc. 

In the LWFA and its variations the laser pulse trav¬ 
elling in the plasma leaves a wakefield of plasma waves 
behind; a bunch of electrons (either externally [TT] or 
self injected [12]) can be accelerated “surfing” one of 
these plasma waves and exit the plasma sample just be¬ 
hind the pulse, in the same direction of propagation of 
the latter (forward expulsion). In Ref. [13] a new LPA 
mechanism, named slingshot effect, has been proposed, 
in which a bunch of electrons is expected to be acceler¬ 
ated and expelled backwards from a low-density plasma 
sample shortly after the impact of a suitable ultra-short 
and ultra-intense laser pulse in the form of a pancake 
normally onto the plasma (see fig. [^. The surface elec¬ 
trons (i.e. plasma electrons in a thin layer just beyond 
the vacuum-plasma interface) first are all displaced for¬ 


ward (with respect to the ions) by the ponderomotive 
force Fp := {—e{'^ x B)^) generated by the pulse, leav¬ 
ing a layer of ions completely depleted of electrons (here 
() is the average over a period of the laser carrier wave, 
E,B are the electric and magnetic fields, v is the elec¬ 
tron velocity, c is the speed of light, z is the direction of 
propagation of the laser pulse); Fp is positive (negative) 
while the modulating amplitude Cg of the pulse respec¬ 
tively grows (decreases). These electrons are then pulled 
back by the longitudinal electric force F^ = —eE^ ex¬ 
erted by the ions and the other electrons, and leave the 
plasma. [In the meanwhile the pulse proceeds deeper 
into the plasma, generating a wakefield.] Tuning the 
electron density in the range where the plasma oscilla¬ 
tion period T^[28] is about twice the pulse duration r, 
we can make these electrons invert their motion when 
they are reached by the maximum of Cg , so that the neg¬ 
ative part of Fp (due to the subsequent decrease of e^) 
adds to Fg in accelerating them backwards; thus the to¬ 
tal work W = Jq dt Fp{v^) done by the ponderomotive 
force is maximal [29]. Provided the laser spot size R is 
sufficiently small a significant part of the expelled elec¬ 
trons will have enough energy to win the attraction by 
ions and escape to infinity. 

Very short r’s and huge nonlinearities make approx¬ 
imation schemes based on Fourier analysis and re¬ 
lated methods (slowly varying amplitude approximation, 
frequency-dependent refractive indices,...) unconvenient. 
On the contrary, in the relevant space-time region a MHD 
description of the impact is self-consistent, simple and 
predictive (collisions are negligible, and recourse to ki¬ 
netic theory is not needed). Here we develop and im¬ 
prove the 2-fluid MHD approach introduced in [T3j [T4) 
and apply it to determine a broad range of conditions 
enabling the effect, as well as detailed quantitative pre¬ 
dictions about it (a brief summary is given in [T511T6]). In 
section [H] we study the plane problem {R = oc) and show 
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1. Free transverse EM wave + 2. Boost by the ponderomotive force exerted 

plasma at rest (plane idealization) by the EM wave (plane wave idealization) 



3. The electric force due to the separation 
of charges boosts the electrons backwards: 
like a SLINGSHOT (plane wave idealization) 


4. Since R < “, the Coulomb attraction by 
ions 0 as Ze ^ , and allows z^ “ 


FIG. 1. Schematic stages of the slingshot effect 


that for sufficiently low density and small times (after the 
impact) we can neglect the radiative corrections [back- 
reaction of the plasma on the electromagnetic (EM) field 
and determine the motion of the surface electrons 
in the bulk by (numerically) solving a single system of 
two coupled first order ordinary differential equations of 
Hamiltonian form, if the initial density is step-shaped, 
or a collection of such systems, otherwise; the role of 
‘time’ is played by the light-like coordinate ^ = ct — z. 
The rough model of m considered only step-shaped no 
and was based on neglecting: during the forward mo¬ 

tion, Fp during the backward motion of the electrons; 
the estimates could be considered reliable only for very 
low, unrealistic Here no needs no longer to be so low, 
nor step-shaped, as in [13], because we take F^^Fp in 
due account during the whole motion of the electrons. In 
section [TTT| we heuristically modify the potential energy 
outside the bulk to account for finite R and determine 
a i?-range such that the motion of the surface electrons 
coming from some inner cylinder be (by 

causality) well approximated by the solution of the cor¬ 
respondingly modified Hamilton equations; we then find 
which electrons indeed escape to infinity and estimate 
in detail their final energy spectrum, collimation, total 


number, charge and energy. To be specific, in section IV 


we specialize predictions to potential experiments at the 
FLAME facility (LNF, Frascati) or the ILIL laboratory 
(INO-CNR, Pisa). We welcome 3D simulations and ex¬ 
periments checking these predictions; the experimental 
conditions are at hand in many laboratories. In section 
|V|we discuss the results, the conditions for their validity 
and draw the conclusions. 


As a context remark, we recall that relatively simple 
2-fiuid magnetohydro dynamic models can be used also to 
describe the complicated physics of the impact of very in¬ 
tense and short laser pulses on overdense solid targets. If 
the density gradient of the target is sufficiently steep, the 
massive displacement of electrons (induced by the pon¬ 
deromotive force) with respect to ions (named snowplow 


in [HI uni) produces a longitudinal electric force which 
may accelerate also protons or other light ions, either 
backward or forward, by the socalled Skin-Layer Pon¬ 
deromotive Aeeeleration m or Relativistieally Indueed 
Transpareney Aeeeleration dHuin] mechanisms. 


The 2-fluid magnetohydrodynamic framework 


The set-up is as follows. We assume that the plasma is 
initially neutral, unmagnetized and at rest with electron 
(and proton) density equal to zero in the region z < 0. 
We describe the plasma as consisting of a static back¬ 
ground fluid of ions (the motion of ions can be neglected 
during the short time interval in which the effect occurs) 
and a fully relativistic collisionless fluid of electrons, with 
the “plasma -1- EM field” system fulfilling the Lorentz- 
Maxwell and the continuity equations. We show a poste¬ 
riori that such a MHD treatment is self-consistent in the 
spacetime region of interest. We denote as cCe(L f^he 
position at time t of the electrons’ fluid element initially 
located at X = (X, V, Z), and for each fixed t as Xe(t, x) 
the inverse map [x = {x^y^z)]. For brevity, we refer: to 
such a fluid element as to the “X electrons”; to the fluid 
elements with arbitrary X, Y and specified Z, or with X 
in a specified region f], respectively as the “Z electrons” 
or the “D electrons”. We denote as m, rig, Ve the electron 
mass, Eulerian density and velocity and often use the di- 
mensionless fields l3e = Ve/c, 

= \/l—/3^ = i/l+Ug. The equations of motion are 


dtXe{t,X) = Ve[t,Xe{t,X)\ 


( 1 ) 


in CGS units {d/dt = dt+VeXx is the electrons’ material 
derivative) and the initial conditions are Pe{0,X) = 0, 
Xe{Q,X)=X for Z > 0. The Lagrangian fields depend 
on t,X, rather than on t,cc, and are distinguished by 
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a tilde, e.g. n^{t^X) = ne\t^Xe{t^X)]. The continuity 
equation dUe/dt^UeVx' Ve = ^ follows from the local 
conservation of the number of electrons, which amounts 
to 


he{t,X)dei[^^=m{X)=he{0,X). ( 2 ) 

We assume that no is independent of X, Y and, as said, 
vanishes if Z < 0; also as a warm-up to more general 
Z-dependence, we start by studying the case that it is 
constant in the region Z>0: h^(Z)=no^(^), where 
0 is the Heaviside step function. We consider a purely 
transverse EM pulse in the form of a pancake with cylin¬ 
drical symmetry around the z-axis, propagating in the 
positive z direction and hitting the plasma surface ^ = 0 
at t = 0. We schematize the pulse as a free plane pulse 
multiplied by a “cutoff” function Xh{p) which is approx¬ 
imately equal to 1 for p = < R and rapidly goes 

to zero for p>R (with some finite radius i?, see fig. Si) 

E\t,x) = Y{ct-z) Xh{p), = zxFY (3) 

[in particular we consider Xaip) = ^{R~p)]'-) fhe ‘pump’ 
e^(^) vanishes outside some finite interval 0<<^</ [30] . 


II. PLANE WAVE IDEALIZATION 


In the plane problem {R = oc) the invertibility of Xe : 
X X for all fixed t amounts to Ze{%Z) being strictly 
increasing with respect to Z for all t. Eq. © becomes 

['^e^]it,Z)=^{Z) ne{t,z)=n^[ZJf.,z)]^{t,z). 

(4) 

Regarding ions as immobile, the Maxwell equations imply 
m that the longitudinal component of the electric field is 
related to N{Z) = Jq dZ' uq^Z') (the number of electrons 
per unit surface in the layer ^<Z' <Z) by 

E^{t,z)=^T:e{N{z)-N[Z,{t,z)]]. (5) 

We partially fix the gauge [14] imposing that the trans¬ 
verse (with respect to z) vector potential itself is in¬ 
dependent of and hence is the physical observ¬ 

able A^(t^z) = —J^dt'cE\t\z); then cE^ = — 

B = B^ =zAdzA^. As known, the transverse component 
of the Lorentz equation implies Pg — ^A^const on 
the trajectory of each electron; this is zero at t = 0, hence 
Pg = mcUg = eA^/c. Hence is determined in terms of 
A^. As in mi, we introduce the positive-definite field 

Se=le-K^ (6) 


which we name electron s-factor. 

covered from Ug , Se through the formulae (44) of [T4| : 


7e=- 




2Se 


2Se 


Up 





FIG. 2. Past (light brown) and future (purple) causal cones 
Dx^T] supports of (light green) and no{z) (anthracite). 


Remarkably, all of Q are rational functions of uf , Sp (no 
square roots appear). Moreover, fast oscillations of uf 
affect but not Sp [see the comments after ([T^]. 

Eor these reasons it is convenient to use uf , Sp instead of 
Ug , i4g as independent unknowns. The evolution equation 
of Sp (difference of the ones of Jp^u^; the former is the 
scalar product of with ppj^pmfc^) reads 

7e^ = —Se + (5t + ca^)tiy. (8) 

dt me 


The Maxwell equation for A^ takes the form {dQ-d‘^)A^-\- 
A^47re^np/mc‘^jp = 0; eq. ^ with R = oo implies 
A^ (^ z) = cP (ct—z) for t < 0, where a^(^) = 

Using the Green function of the D’Alembertian 9 q— 9^, 
abbreviating x = (t,^), these equations can be equiva¬ 
lently reformulated as the integral equation (42) of [14] 


A^{t,z) — of{ct — z) = — f dt' dz' 

J DxCiT 


2TTe^np 

mc7e 


(x') (9) 


Dx = {(t'^ z') 11' <t, \z — z'\ <ct — ct'}^ T = {x I <ct} 

The past, future causal cones Dx^T, the supports of 
A^,no{z), and their intersections are shown in fig. 

Eor t < 0 DxPiT is empty, and the right-hand side of Hi 
is zero, as it must be. Below we shall analyze the con¬ 
sequences of neglecting it also for small t, and determine 
the range of validity of such an approximation. 


II.1. Motion of the electrons 


Let u^{^) = ea^{^)/mc^, 

F|(2,Z) = -47re2 |V(0)-iV(Z)|. 

Ff{t,Z) = Ff[zp{t^Z)^Z] is the longitudinal electric force 
acting on the Z electrons at time t; it is conservative. 
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as it depends on t only through Ze{t^Z). The approxi¬ 
mation {^z) = cx^ {ct — z) implies {% z) = {ct — z ), 
and the last term of (|8) vanishes. Replacing <§ in the 
Lagrangian version of (^, we find for each Z > 0 the 
equation jedtSe = —s^F^jmc. The initial condition is 
Se(0, Z) = 1. The other equation to be solved is Q 2 with 
the initial condition Xe{0^X) = X. By one is thus 
led to the Cauchy problems (parametrized by Z>0) 


di 


Ze-Z _ lFv[ct-Ze{t,Z)]-sl 


lFv[ct-Ze{t,Z)]Fsl ’ 


dtSe = (11) 

7emc 


^e(0, Z) — Z — 0, 


SeiO,Z) = l. 


( 12 ) 


Xg(t,X) is obtained from the solutions of ( 11||12 ) using 

0, 0: 

x^{t,X) = X^ + [ dt'c/3^[t',Ze{t',Z)]. (13) 

Jo 


For all fixed Z the map t C{t^Z) = ct — Ze{jZ) is in¬ 
vertible, because the speed of electrons is always smaller 
than c. We can simplify 0 by the change of variables 
(t,Z) 1 -^ making the argument of v an indepen¬ 

dent variable. Denoting the dependence on (^, Z) by a 
caret [e.g. s{^,Z) = Se{t,Z)] and introducing the displace¬ 
ment from the initial position A(^,Z) = Ze(^,^) —we 
find = (7e/cse)9t, and 0 becomes 

(the prime means differentiation with respect to ^). For 
^<0 'z;(^) = 0, A, 5 remain constant, and we can replace 
the initial conditions A(—Z,Z)=0, s(—Z,Z) = 1 by 


A(0,Z)=0, 


s(0,Z) = l. 


(15) 


An alternative derivation of ( |l4p5 ) with a deeper insight 
on the role of the s-factor is given in [20]. In the zero 
density limit N{Z) = 0, s = i, diiliit is integrable, and 
all unknowns are determined explicitly from [nmi]). 
As > 0, even if u\ v oscillate fast with integrating 
(14) makes relative oscillations of A much smaller than 


those of V and those of s much smaller than the former; 
hence, s is practically smooth, see e.g. fig. Setting 
g = A, p= —s, for each fixed Z (14) are the Hamilton 
equations (with ‘time’ q' = dHjdp, p' = —dH/dq of 
a system with Hamiltonian H{q,p^^; Z) = H (g, —p, Z), 


H{A,s,^;Z) ^ 7 ( 5 , 0 + ^(A;Z), 7 ( 5 , 0 ^ 

U{A-,Z) = ^ \M{Z+A)-M{Z)-N{Z)A 


N{Z) = J^dZ'N{Z') = J^dZ'm{Z') {Z-Z'). 
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(16) 


Defining U we have fixed the free additive constant so 
that Z//(0,Z) = 0 for each Z; is positive definite. 

Below we shall abbreviate P(^; Z) = (A(^; Z), Z )). 


The right-hand side of (14)2 is an increasing function 
of A, because so is N{Z). As v{^) is zero for ^ < 0 
and positive for small ^ > 0, then so are also A(^,Z) 
and 5(^,Z) —1. Both keep increasing until A reaches a 
positive maximum A(^, Z) at the ^ = (^(Z) >0 such that 


A'(C,Z)=0 s\^,Z) = l+viO (17) 


(note that ^ < I if v{l) = 0). We shall denote as 
C = A[^(0),0] the maximum penetration of the Z = 0 


electrons. For ^>^(Z) A starts decreasing; 5 reaches a 
maximum at the ^ = ^r(Z) such that A(^^,Z)=0 (i.e. at 
^ = ^^(Z) the Z electrons have regained their initial z). 
Both decrease for ^ > ^r(Z), until 5 becomes so small, 
and the right-hand side of (14)i so large, that first A, 
and then 5 — 1, are forced to abruptly grow again to pos¬ 
itive values. This prevents 5 to vanish anywhere, consis¬ 
tently with In ^-intervals where v{^)=Vc=coiist^ H 
is conserved, and all trajectories P(^; Z) in phase space 
(paths) are level curves H{A,s;Z) = h{Z)^ above the 
line 5 = 0, integrable by quadrature [22]. For Z = 0 the 


paths are unbounded with A(^, 0) ^ — 00 as ^ ^ oc. For 
Z > 0 the paths are cycles around the only critical point 
C = (A, 5 ) = (0, y/lFVc) (a center); therefore for ^ > / 
v{^)=v{l), and these solutions are periodic. There exists 
a Z 5 > 0 such that: the paths Z) with Z < Zij cross 
the A = —Z line twice, i.e. go out of the bulk and then 
come back into it; the path P(^; Z 5 ) is tangent to this 


line in the point {A, s) = {—Zb, ^/lFv(J)) (where A' = 0); 
the paths P(^; Z) with Z> Z}j do not cross this line. For 
Z < Zi) let ^ex(Z) be the first positive solution of the 
equation A(^,Z) = —Z, i.e. at ^ = ^ex(^) the Z electrons 
exit the bulk: 


Zel^ex(Z),Z]=0. (18) 


The function ^ex(Z) is strictly increasing if dzZe>0. 
For any family P{^]Z) of solutions of ( T4p5 ) let 


u = - 


~Ys ^ 


7=- 


~Ys ^ 




.(e,X) = X + Y{^,Z), Y{^,Z)=ldy (19) 


KvZ) 


cf{^, Z)=Z+'={^,Z), E{^,Z) = jdy =C+A(^,Z) 

Jo s{y,Z) 


(note that = A). The so defined 'U,7,Xe are the so¬ 
lutions - expressed as functions of X - of all equations 
and initial conditions m- Note that Xe,t can be ob¬ 
tained also solving the system of functional equations 


^ = ct—z^ E{ct—z^Z) = ct—Z, X —X = Y{ct—z^Z) 

( 20 ) 

[by ( [T^ the second is actually equivalent to the 2 :- 
component of the third] with respect to t^x. Clearly 
S(^,Z) is strictly increasing and invertible with respect 
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to ^ for all fixed Z. Solving ( [2Q| ) with respect to x 
(resp. as functions of t, X (resp. of t^x) and re¬ 

placing the results in ^,7, s,... one obtains the solutions 
in the Lagrangian (resp. Eulerian) description: in par¬ 
ticular one finds (generalizing [14]) 


|(t, Z)=E-\ct-Z,Z), Xe{t,X)=X^Y |(t, Z),Z 


Zq (i^Z^ — Z “h X 


i{t,z),z =ct-i{t,z), 

Ue{t,Z)=u\^{t,Z),Z 


Se{t^Z) = S ^{t^Z)^Z 
X^{t^x) = X^ —Y^[ct — Z^Ze{t^Z 
Ue{t^z)=u[ct — Z^Ze{t^z)\. 


( 21 ) 


Indeed, it is straightforward to check that 
(^Ze{t, Z),Se{t, Z)) is the solution of ( [ll][T^ and 
Pe(t^x)=mcUe(t^z)^ Xe(t^X) of the PDE’s(0 With the 
initial conditions pe{0, X) = 0, Xe{0, X) = X for Z>0. 

Erom ( [T^ , ([1^5, the times of maximal penetra¬ 
tion and of expulsion of the Z electrons are 


= = ( 22 ) 


Deriving (21) and the identity y = 
obtain a few useful relations, e.g. 


7-1 


{y,z)Z 


we 


^ —s dA 


dZe 

dZ ~ j dZ 


’ dZ 


-\v,Z) 

dZ, _ 

7 


dz sdzZe 



^=S-i(ct-Z,Z) 


(23) 


By (23), 9^Ze = H-9^A>0 is thus a necessary and suffi¬ 


cient condition for the invertibility of the maps Ze\Z\-^ z^ 
Xe'.X\-^x (at fixed t), justifying the hydrodynamic de¬ 
scription of the plasma adopted so far and the presence 
of the inverse function Z^it^z) in (21). Einally, from 
(23) we find also 


ne{t,z)=no[Ze{t,z)] 


7 


sdzZe 


i^,Z)={ct-Z,Zeit,z)) 


(24) 


We can test the range of validity of the approximation 
A^{^z) = cY{ct — z) by showing that the latter makes 
the modulus of the right-hand side of (§ much smaller 
than cY{ct — z) on D = {{t,z) | 0 < ct —2: < 0 < 

ct + z <^ex{ZM)} {Zm is defined below), or equivalently 
[multiplying by elmc? and using 


for x = {t,z)e D z) \ <C |i/(ct —2))|, (25) 

■'anr 


actually, it suffices to check this inequality on the world¬ 
lines of the expelled electrons. 


II.2. Auxiliary problem: constant initial density 


As a simplest illustration of the approach, and for 
later application to a step-shaped initial density, we first 
consider the case that no{Z) = no. Then is 
the force of a harmonic oscillator (with equilibrium at 
Zg = Z) F^{ze^Z) = —ATTnoe^[ze — Z] = —47rnoe^A; the 


Z-dependence disappears completely in ( 14|[T5 ), which 
reduces to the auxiliary Cauchy problem 


A(0)=0, 5(0) = 1, (26) 


where M = Aire^no/mc^. The potential energy in ( [T^ 
takes the form U{X^Z) = MA^/2. Problem (26), and 
hence also its solution (A(^),5(^)), the value of the en¬ 


ergy as a function of ^ and the functions defined in (19), 


are Z-independent. It follows S^A = 0 and by (23) the 


automatic invertibility of Ze{t,Z); moreover, the inverse 
function Ze{t,z) has the closed form 

Ze{t,z) = ct—E{ct — z) = z—A{ct — z) (27) 

[here 5(^) =^+A(^)], what makes the solutions ( [^ of the 
system of functional equations (20), as well as those of 


(|y , completely explicit in terms of E and the inverse 
only. As a consequence, all Eulerian fields depend on t, z 
only through ct—z (i.e. evolve as travelling-waves). In fig. 
I^left we plot some solution of ( [2^ . If 'z;^) ='r’g = const 
all paths P{^;Z) are cycles around C (fig. [fright), corre¬ 
sponding to periodic solutions. Within the bulk electron 
trajectories for slowly modulated laser pulse like the ones 
considered in section [IV| are tipically as plotted in fig. 
in average they have no transverse drift, but a longitu¬ 
dinal forward/backward one. Eig. shows a couple of 
corresponding charge density plots. 


III. 3-DIMENSIONAL EFFECTS 

We now discuss the effects of the finiteness of R. For 
brevity, for any nonnegative r, L we shall denote as Cr 
the infinite cylinder of equation p<r, as Cf: the cylinder 
of equations p<r, 0<z<L. The ponderomotive force of 
the pulse will boost forward (as in fig. only the small- 
Z electrons located within (or nearby) Cr. These For¬ 
ward Boosted Electrons (FBE) will be thus completely 
expelled out of a cylinder which will reach its maximal 
extension Cr around the time f(0) of maximal longitudi¬ 
nal penetration C = A[f(0),0] of the Z = 0 electrons. The 
displaced charges modify E. By causality (see appendix 
[a] ), for X near the z axis E{t^x) is the same as in the 
plane wave case for t<t( 0)+i?/c, and smaller afterwards. 
We choose no, so that they fulfill 

r = i?-fc^0(rfe.-O>O (28) 

and condition ( [^ for all x = (t, x) such that t < t(0)+i?/c; 
here t = t(0), tex='f^ex{^) a-re the times of maximal pen¬ 
etration and of expulsion from the bulk of the Z = 0 
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FIG. 3. (a) Solution of (26) for Ml^ = 26 and the as 
in section IV of average intensity /= lO^^W/cm^. (b) Paths 
P(^; Z) around the center C for = 26, Vc^O. 


electrons [see ([2^ ]. [As grows from zero the right- 
hand side of (|25p does as well, whereas t(Q), tea;(Q) — t(Q) 
decrease]. In appendixj^we show that conditions ( [^ re¬ 
spectively ensure that these FEE, at least within an inner 
cylinder p<r < R: i) move approximately as in section 
\n\ until their expulsion; ii) are expelled before Lateral 
Electrons (LE), which are initially located outside the 
surface of and are attracted towards the ^-axis (see 
fig. ^3), obstruct their way out. Eor the validity of our 
model we must a posteriori check also that the expelled 
Cr electrons remain in 

i.e. their transverse oscillations Axe are <C R. (29) 

In the plane model the expelled Z > 0 electrons can¬ 
not escape to 2; ^ — oc because are decelerated by the 
constant electric force Z) = 47 Te‘^N{Z) >0, see (10). 
The real electric force > 0 acting on the Cr electrons 
after expulsion is generated by charges localized in Cr, 
hence F^^ as ^ — 00, and the escape of expelled 

electrons is no more excluded a priori. Moreover, since 
F^{t,0) = 0, it should be also F^^{t,0) = 0, allowing the 
escape of the Z = 0 electrons; by continuity there will 




I 

-0 

I 

I 


f 


FIG. 4. Normalized charge density plot l — rie/no under the 
same conditions as in fig. [Rafter about 25 fs (up) and 70 fs 
(down), i.e. resp. before and after the maximal penetration 
time t(0) = 51fs; in the latter picture the electrons travelling 
backwards make light yellow-blue striped the region between 
the yellow and the white-blue striped ones. 


exist some positive Zm < Zi) such that the C^^ electrons 
escape to infinity. We stick to estimate F^^ on the ^-axis 
electrons; we assume that after the pulse has overcome 
them, they move along the z-axis. Actually this will be 
justified below if u^{l) 0, which in turn holds if, as 

usual, I ^ A [see the comments after (33)]. In fig. a) 
we schematically depict the charge distribution expected 
shortly after the expulsion. The light blue area is occu¬ 
pied only by the X G C^^ electrons. The orange area 
is positively charged due to an excess of ions. For any 
Z-electrons moving along the ^-axis consider the surfaces 
Sq, 81,82 occupied at time t by the X' G Cr electrons re¬ 
spectively having Z' = 0, Z, Z 2 {Z), where Z 2 {Z) is defined 
by the condition N{Z 2 ) = 2 N{Z), which ensures that the 
electron charges contained between 80 , 81 and 81,82 are 
equal (in the figure 80 , 81,82 are respectively represented 
by the left border of the blue area, the dashed line and 
the solid line). The longitudinal electric force F^'^ acting 
at time t on this Z-electron is nonnegative and can be 
decomposed and bound as follows m- 

0 < Ff{t,Z) = -eEi{t,Z)-eEi{t,Z) < F:,[A{t,Z),Z]. 

Here E^{t,Z) stands for the part of the longitudinal 
electric field generated by the electrons between 80 , 82 ; 
since those between 80 , have by construction the same 
charge as those between 81 , 82 ^ but are more dispersed, 
it will be —eE^_(t, Z) < 0. The part —eE^{t,Z) of F^^ 
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FIG. 5. a) schematic picture of the expected charge distri¬ 
bution shortly after the expulsion (long arrows) of surface 
electrons; short arrows represent the inward motion of the 
lateral electrons; b) simplihed charge distribution generating 
the effective potential energy. 


generated by the ions and the remaining electrons (at the 
right of S 2 ) will be smaller than the force generated 
by the charge distribution of fig. [^b), where the remain¬ 
ing electrons are located farther from (0,0, ^e) (in their 
initial positions X', not in the ones at t) and hence gen¬ 
erate a smaller repulsive force. This explains the second 
inequality in the equation. In appendix we show that 
for Ze = Z-\-A<0 


F!r{^,Z) 

27re^ 


= 2 N{Z) - [ 
Jo 




dZ^ 


0 y(Z'-^e)2+r2' 


(30) 


Commendably, is conservative, nonnegative and goes 
to zero as A ^ — 00 , while it reduces to zero for Z = 0 and 
to A7Te^N{Z) as r^oc, as F^ in ( 10 > 3 ; it becomes a func¬ 
tion of t (resp. through A(t, Z) [resp. A(^,Z)] only. 
We therefore modify the dynamics outside the bulk re¬ 
placing F^ by , or equivalently U hy Up in ( [l^ , where 
Ur is continuous and equals U for 2 (e = ^+A>0, and the 
potential energy (B2) associated to F^^ for Ze = Z-\-A<0] 
there Ur is a de creasing function of A with finite left 
asymptotes ( |B3[ ). We will thus underestimate the final 
energy of the electrons, because F^^ is larger than the 
real electric force F^^ decelerating the electrons outside 
the bulk; this makes our estimates safer. In fig. [^we plot 
suitably rescaled U and Ur for no{Z) = noO{Z). After 


the pulse is passed we can compute 7 as a function of 
A, Z using energy conservation mc^jFUr{A,Z) =const. 
For the expelled electrons the final relativistic factor 
7 j(Z) = 7 e(A = —oc,Z) is the decreasing function (B4). 
The maximum of yf{Z) is 7 m =7/(0)- Let Zm < ^6 be the 
Z fulfilling jf{Z) = 1. The estimated total number Ae, 
electric charge (in absolute value) Q, and kinetic energy 
E of the X G escaped electrons are thus 


pZm ( 31 ) 

E^mTTC^r^/ dZno{Z)[yf{Z) — l]. 

Jo 

The number of escaped X' G electrons with Z < 

Z' < ZFdZ is estimated as 7 rr‘^no{Z)dZ, that with rel¬ 
ativistic factor between 7 and 7 + ^7 is estimated as 
dN = 7rr‘^[rio{Z)/\dyf/dZW^^^^^^^dy^ where Z{y) is the 
inverse of 7 / (Z) (a strictly decreasing function, see ap¬ 
pendix]^. Hence the fraction of escaped electrons with 
final relativistic factor between 7 and 7 + 6/7 is estimated 
as 1 /( 7 ) 6 / 7 , where 


2 /( 7 )= 


j_^ 


1 no{Z) 
N{Zm) \<ilf/dZ\ 


Z=Z(7) 


(32) 


determines the associated energy spectrum. As (a^(0 = 
by 0 the final transverse deviation of the 
escaped electrons will be 


u3 u3 


(33) 


where uj = u\l). This is an increasing function of Z, 
because 7 /(Z) is decreasing. If X<t^l then uj (see 
next section), and (33) is negligible unless Z Zm • 


III.l. Step-shaped initial density 


If no{Z)=no0{Z) then N{Z)=no6{Z)Z, and for Z>0 


F^^{A,Z) f- 2 A (elastic force), 


Zf=^0 


(34) 


27rnoe2 [2Z+ ^Z+A^+r^ - ^Z-A^+r^, Ze<0: 

Since the first expression is the same as in the case 
rio{Z) = n o, th e motion of the Z-electron will be as in 
subsection II.2 until ^ = ^ex(Z). The second expression 
goes to the constant force dyrnoe^Z as r ^ oc, as ex¬ 
pected. The motion for ^>^ex{Z) will be studied in de¬ 
tail in [ 22 ] ; we plot the graphs of a typical solution (until 
the expulsion) in fig. and a few corresponding electron 
trajectories in Fig. We can readily understand that it 
will be 9^Ze(^,Z)> 0 for all ^ and 0 < Z < Zm, since this 
holds for ^< ^ex{Z) [by the comments following ( [^ ] , and 
both ^ex{Z) and the decelerating force Fg^(A,Z)^utside 
the bulk) increase with Z, while the speed of exit from 
the bulk decreases with Z, whence the distance between 





























FIG. 6. Rescaled longitudinal electric potential energies u = Uj^ Ur = Ur j for (left) idealized plane wave 
Rjl^oo or (right) for R// = 0.85, plotted as functions of A for ZjZM — 0, .2, .4, .6, .8,1; the horizontal dashed lines are the left 
asymptotes of Ur for the same values of ZjZM- Here the initial electron density is step-shaped: no{Z) = noO{Z). 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 



electrons with different Z increases with t. The Z 5 in¬ 
troduced before ( p^ is now the solution of the equation 
MZ^/2 — i.e. the Z corresponding to the 
zero longitudinal velocity and the final value of the energy 
h after the interaction of the pulse; one can determine h 
evaluating at^ = /, /z= |{s(/)+[l-K;(/)]/ 5 (/)+M[A(/)]^}. 

Hence, 


-1 2 


Zb = \[A{l)]^+ sil)-^l+v{l) /2Ms(l). (35) 


7 / {Z), v{l ) admit rather explicit forms ( |B7[ ), ( |B 8 [ ). In 
section |IV| we plot spectra 1 /( 7 ) corresponding to sev¬ 
eral no and intensities. Moreover, Q = Trr^eno^M, 
E = 'Kr‘^nQmc^j^^dZ[^f{Z) — l]. Finally, if Cetc(0)<^ then 
in ( [^ becomes 


ct—z ct-\-z 







(36) 


IV. NUMERICAL RESULTS 



FIG. 7. (a) Laser pulse of a verage intensity I = lO^^W/cm^ 
and shape as in section EL with I = 18.75/im. (b-c) Gor- 

responding solution of ( 14] 15 ) for initial density n^{Z) = 
no0{Z), with no = 21 x 10^cm“^ (i.e. M/^ = 26). 


We assume for simplicity that the pulse is a slowly 
modulated sinusoidal function linearly polarized in the x 
direction: e^(^) = es{^)xcosk^, the modulating am¬ 
plitude €s{^) > 0 is nonzero only for 0 < ^ < /, and 
slowly varies on the scale of the period A = 27r/k I, 
i.e. Aje'^l <C | 6 s| on the support of e^. Integrating by parts 
we find = xes{^)lsmk^)/k^O{l/k‘^) [ 20 ] and, in 

terms of the rescaled amplitude w{^) = ees{^)/kmc^, 

sm{k^), v(^)~«;2(C) sin^(fc^), (37) 

where ac^b means a = b+0{l/k‘^). Note that, as es(^) = 
0 for ^>/, this implies uj = u\l)c^0, as anticipated. 

If we approximate as Xr{p) = ^{R~p) the cutoff function 
in ® , the average pulse intensity on its support is I = 
cSpKR^^l. Here £ is the EM energy carried by the pulse. 





d 2 rl d 2 rl 

( 38 ) 
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FIG. 8. Trajectories gone in ca. 150 fs by electrons initially located at Z/Zm — 0, 0.25, 0.5, 0.75,1 under conditions as in fig. 


High power lasers produce pulses where A^l/iui and 
is approximately gaussian, es(^) (xexp[—(^ —^o)^/2cr]; a 
is related to the fwhm (full width at half maximum) 
of by cr = /'^/41n2. If initially matter is composed of 
atoms then es{ct — z) can be considered zero where it is 
under the ionization threshold, because the pulse has not 
converted matter into a plasma yet. Hence we adopt as 
a modulating amplitude 6s(^) the cut-off Gaussian 


eg (0 = ^3 exp 


-g-V2)^ 


2cr 


/2 


,2 16 Vin 2 e 
bi = -; 


R‘^1' 


R = 


Oimi-a ( 39 ) 


\/In2 


In 


41n2^ 
\/ln 2 mc^ E (eA)^ 


Ui ^/nl' {'kRtyic^Y 


where Ui is the first ionization potential (for Helium Ui ^ 
24eH); the formula for 6^ follows replacing the Ansatz 
([^1 in [neglecting the tails left out by the cutoff 
^(^^(/ —Numerical computations are easier if we 
adopt m as €s{^) the following cut-off polynomial: 


bn r 


2 


ep{C)^f[i-{2aip-ir\ e{O0{ip-O, (40) 

bpjp are determined by the requirement to lead to the 
same fwhm and £: bp = 5040 £/Ip and lp = bl'I2. 

We now present the results of extensive numerical sim¬ 
ulations based on the experimental parameters available 
already now at the FLAME facility [23] or in the near 
future at the ILIL facility [32]: I' 7.5/im (implying 
/p = 18.75/im), (implying kip = 27Tlp/Xc^ 147), 

f = 5J, and R tunable by focalization in the range 
10“^ ^ 1 cm. We model the electron density: first as the 
step-shaped one ri^{Z) = no0{Z) (this allows analytical 
derivation of more results); then as a function smoothly 
increasing from zero to the asymptotic value no, with 
substantial variation in the interval d<Z<L = 20/im (as 
motivated by experiments, see section]^, more precisely 
no{Z) =no0{Z)tdinh.{Z/L). We have numerically solved 
the corresponding systems ( T4p5 ) and proceeded as in 
section |HI[ for R = 16,15, 8,4, 2,1/im [resp. leading to 
average intensities //10^^(W/cm^) 1,1.1,4,16, 64, 255], 


no in the range lO^^cm ^ < no < 3 x lO^^cm ^ and 
Z<Zm] all results follow from these solutions. 

In fig. [^left we plot the maximal final relativistic 
factor 7 m of the expelled electrons as a function of no, 
with the above values of I and no(^) = '^o^(^); each 
graph stops where no becomes too large for conditions 
(28)i, or (29) to be fulfilled and is red where con¬ 
dition (28)2 is no more fulfilled. The latter prevents 


collisions with the LE and becomes superfluous if the 
target is a solid cylinder of radius R (since then there 
are no LE) [33]; the / = 64, 255 x lO^^W/cm^ graphs 
are plot green for densities corresponding to the light¬ 
est solids (aerogels) available today. As expected [13]: 1) 
as no ^ 0 7 m — 1 oc 2) each graph 7 m(^ o; I) has 

a unique maximum = ^Mi^OMp) at noM ^ f^o, 

where no is the density making ^(0) = //2, namely such 
that the Z = 0 electrons reach the maximal penetration 
C = A(//2,0) when they are reached by the pulse max¬ 
imum. The dependence of 7 m on no is anyway rather 
slow. The striking 7 mm (^) cx: / behaviour shown in fig. 
[^ up-center hints at scaling laws and will be discussed 
elsewhere. In figures we plot sample spectra ^(7) 
for //10^^(W/cm^) 1,4,16, 64 and tlq compatible with 

( [2^ , ( [2^ , ( [2^ . In table 1^ we report our main predictions 
for the same / (equivalently, R) and tiq. The final ener¬ 
gies of the expelled electrons range from few to about 15 
MeV. The spectra (energy distributions) are rather flat 
for the step-shaped densities, albeit they become more 
peaked near 7 m as no grows; if no{Z) grows smoothly 
from zero to about the asymptotic value no in the inter¬ 
val 0 < Z < L ^ 20/im, they can be made much better 
(almost monochromatic) by tuning L. The collimation 
of the expelled electron bunch is very good, by (33); in 
all cases considered in table |I| we find deviations Pf/Pf 
of 1 d- 2 milliradiants for the (p, Z) = (0, 0) and 4d-10 mil- 
liradiants for the (p, Z) = (0 , 0.9Zm) expelled electrons. 

We now discuss the conditions guaranteeing the valid¬ 
ity of our model. The comments after (34) show for all ^ 


the invertibility of the maps -): Z z in the inter¬ 
val 0<Z < Zm, and therefore the self-consistency of this 
2-fluid MHD model, in the step-shaped density case; nu- 
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FIG. 9. Left: relativistic factor 7 m of the Z — {) expelled electrons (the maximal one) as a function of the step-shaped initial 
electron density no, for few values of the intensity /; the maximum of each graph is denoted as 7mm. Center: 7mm vs. I. Right: 
iF & its correction (5iF along the X=0 electrons’ worldlines for no = 24x lO^^cm”^, 7 = 255 x 10^^W/cm^: is negligible. 






FIG. 10. Sample spectra of the expelled electrons for pulse amplitudes of the form (40 ) with continuous initial electron densities 
^{Z)=ri{)0{Z)ld,-nh{Z/L)^ L = 20/im (graphs a-d), or step-shaped initial electron densities ^ = nQ0{Z) (graphs e-f). The 
values of no and of the average pulse intensity I are the same as in Table 


merical study of the map : Z^z shows that this 

holds true also in the continuous density case. Numer¬ 
ical computations show that (25) is fulfilled at least on 
the Z < Zm electrons’ worldlines, even with the highest 
densities considered here (see e.g. fig. [fright). Finally, 
the data in table show that (28), ( [^ are fulfilled. 


If we choose es{C) as the cut-off gaussian, instead of 
the cut-off polynomial, convergence of numerical compu¬ 
tations is slower, but the outcomes do not differ signifi¬ 
cantly. Sample computations show that choices of other 
continuous ^{Z) lead to similar results, provided the 
function no(Z) is increasing and significantly approaches 
the asymptotic value no in the interval 0 < Z < L ^ 20/im. 


V. DISCUSSION, FINAL REMARKS, 
CONCLUSIONS 


These results show that indeed the slingshot effect is 
a promising acceleration mechanism of electrons, in that 
it extracts from the targets highly collimated bunches of 
electrons with spectra which can be made peaked around 
the maximum energies by adjusting 7?, no; with laser 
pulses of a few joules and duration of few tens of fem¬ 
toseconds (as available today in many laboratories) we 
find that the latter range up to about ten MeV (it would 
increase with more energetic pulses). The spectra (dis¬ 
tributions of electrons as functions of the final relativis¬ 
tic factor 7j), their dependence on the electron density 
and pulse intensity, the collimation and the backward 
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pulse energy wavelength A 

~0.8/im, fwhm 1' 

7.5/im, spot radius Rc^ 

1 ^ 

16/im 


P 

PoP14 

P 

PoP14 

P 

g 

cp 

eg 

cp 

eg 

cp 

eg 

cp 

eg 

pulse spot radius R (ym) 

15 

15 

16 

16 

2 

2 

16 

16 

16 

16 

8 

8 

4 

4 

mean intensity I (lO^^W/cm^) 

1.1 

1.1 

1 

1 

64 

64 

1 

1 

1 

1 

4 

4 

16 

16 

initial el. density no(10^®cm“^) 

.64 

.64 

.64 

.64 

64 

64 

3.2 

3.2 

8 

8 

32 

32 

160 

160 

ratio [tex{0)—t{0)]c/R 

0.8 


.75 


1.3 

1.9 

0.8 

0.8 

0.6 

0.7 

1.1 

1.1 

1.7 

2.1 

ratio r/R 

0.6 


0.7 


1 

1 

0.6 

0.7 

0.8 

0.9 

0.8 

1 

1 

1 

ratio Ax^jR 

.02 


.02 


.14 

.09 

.02 

.19 

.02 

.17 

.05 

.04 

.1 

.06 

maximal relativistic factor 7 m 

2.5 

1.83 

2.3 

1.65 

16 

14 

2.5 

2.4 

2.7 

2.6 

6.6 

5.6 

11 

8.1 

max. expulsion energy iL(MeV) 

1.3 

0.94 

1.2 

0.9 

8.1 

7.2 

1.3 

1.2 

1.4 

1.3 

3.4 

2.9 

5.4 

4.2 

tot expelled charge Q (10“^°C) 

1.7 

3.8 

2.2 

3.2 

3.7 

3.1 

1.9 

2.2 

3.7 

2.9 

3.5 

4 

3.8 

3.5 

tot. exp. kin. energy E(10“^J) 

0.7 


0.7 


16 

12 

0.8 

0.8 

1.6 

1.1 

4.5 

4.0 

8.6 

5.8 


TABLE L Sample inputs and outputs for possible experiments. In the ‘p,g’ columns the initial electron densities are step¬ 
shaped, riZ{Z) = no 0{Z), and the amplitudes are resp. of the gaussian, polynomial forms (39), (40); in the ‘PoP14’ columns 
we report results computed in m with poorer approximation. In the ‘cg,cp’ columns the initial electron densities are the 
continuous ones n^{Z) =no 0{Z) i3inh{Z/L) with L = 20/im, and the amplitudes are resp. of the forms (39), (40). 


direction of expulsion in principle allow to discriminate 
the slingshot effect from the LWFA or other acceleration 
mechanisms. In table [T| and fig. [^we have reported de¬ 
tailed quantitative predictions of the main features of the 
effect for some possible choices of parameters in experi¬ 
ments at the present FLAME, the future upgraded ILIL 
facilities, or similar laboratories. Low density gases or 
aerogels (the lightest solids available today) are targets 
with appropriate electron densities. 

The steepest 2;-oriented density gradient of a gas sam¬ 
ple isolated in vacuum is attained just outside a nozzle 
expelling a supersonic gas jet in the xy plane; across the 
lateral border of the jet the density may vary from about 
zero to almost the asymptotic value no in about L ^ 20/im 
[23] . Hence if we choose a supersonic helium jet as the 
laser pulse target the initial electron density is reasonably 
approximated by the choice ^{Z) = no 6>(Z) tanh(Z/I/), 
and the predictions of table fig. (a-d) are reliable. 
By the way, the values of no considered in table [T] are 
considerably higher than in typical LWFA experiments. 

Step-shaped no(^) are unrealistic approximations of 
densities of gas samples, but reasonable ones of solids 
(for which AZ^A), provided no exceeds 48xl0^^cm“^, 
which is the electron density of aerographene (the light¬ 
est aerogel so far: mass density=0.00016 g/cm^). Sil¬ 
ica areogels, with a wide range of densities from 0.7 to 
O.OOlg/cm^, electron densities of the order of 10^^/cm“^ 
and porosity from 50 nm down to 2 nm in diameter (i.e. 
much smaller than A) have been produced and exten¬ 
sively studied [SUED- Therefore the results of the last 
two ‘p,g’ columns of table I] [and the corresponding spec¬ 
tra, figure [Tq| (e)] are applicable to aerogels, while those 
of the first four are presently only of academic interest. 

The quantitative predictions of our model are based 
on a rather rigorous plane-wave, 2-fiuid magnetohydro¬ 
dynamic model m and simple, but heuristic approxima¬ 


tions for the 3D corrections, which certainly affect their 
liability. We welcome numerical 3D simulations (particle- 
in-cell ones, etc.) to improve the latter. Experimental 
tests are easily feasible with the equipments presently 
available in many laboratories. We welcome experiments 
testing the effect. 
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Appendix A: Finite R conditions 

As known, for any spacetime region V its future 
Cauchy development D^(fD) is defined as the set of all 
points X for which every past-directed causal (i.e. non¬ 
spacelike) line through x intersects V (see fig. [IT] left). 
Gausality implies: If two solutions of the system of dy¬ 
namic equations coincide in some open spacetime region 
V, then they coincide also in D^ifD). Therefore, knowl¬ 
edge of one solution determines also the other (which we 
will distinguish by adding a prime to all fields) in D^(fD). 

In the problem at hand the solutions are exactly known 
for t < 0, i.e. before the laser-plasma interaction begins. 
We use causality adopting: 1. as P a region (see 
fig. right) of equations — e < t < 0 and either p < R 
or z > 0, with some e > 0 (we can take also e = 0 if we 
assign on also the time derivatives of the A^,u); 2 . 
as the known solution the plane one induced (section [n| 
by the plane transverse electromagnetic potential, which 
can be approximated as {t, z) = oA {ct — z) under the 
assumption ( [^ ; 3 . as the unknown solution the “real” 
one induced by the “real” laser pulse A^{t,x), which we 
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FIG. 11. Left: future Cauchy development of a generic domain V. Right: P? (light blue) and (shaded region 

between the blue and light blue hyper surf aces) in {p,z,ct) coordinates (we have dropped the inessential angle (p); worldlines of 
the X = 0 electrons (red) and of a couple of off-z-axis electrons (yellow); the former remain in Il^(Pi) longer. 


approximate as a potential leading to Q. It is easy to 
show that is the union of three regions, resp. of 

equations: a. z>ct>0;h. ct>z>0 and — z‘^ < 

R; c. t>0, z<0 and p^ct<R (see fig. [^. In 
the two solutions coincide, in particular a “real” electron 
worldline remains equal to the plane solution 

worldline Xe{t^X) as long as Xe{t^ X) (V^). 

By continuity, we expect that the two solutions remain 
close to each other also in a neighbourhood of 
This is confirmed by estimates [13] involving the retarded 
electromagnetic potential (in the Lorentz gauge d' A = ^) 


the most dangerous LE (the Z = 0 ones) after the pulse 
has completely entered the bulk is less than R] stands 
for the average p-component of the velocity of these LE. 
By geometric reasons <v^= average ^-component of 
the X = 0 electrons velocity in their backward trip within 
the bulk; our rough estimate v^c^v^l2 = (l{tex—t)2 gives 
(28)2- Eq. (28) is thus explained. 


Appendix B: Finite R energies 


A^(t,x) = Aj (t,x) + Jd 


,3 J^[tr{t,x-x'),x'] 

^ \x — x'\ 


(Al) 


i.e. the general solution of the Maxwell equation \3A^ = 
Airj^ with a current j^(t^x) vanishing for t < 0; here 
tr{t,x — x')=t—\x — x'\/c, A^{t,x) fulfills □Aj=0 (de¬ 
termining the t^—oc behaviour), and E = ^dfA—VA^^ 
B = \/ X A Since the formation of Cr is c omp leted at 
t = t(0), and the ‘information’ [encoded in (Al)] about 
the finite radius of Cr takes a time R/c to go from the 
lateral surface p = R to the Z-axis, then if eq. ([28|)i is ful¬ 


filled the X = 0 electrons (red worldline in fig. 11) move 
approximately as in section [ll| until the expulsion. Simi¬ 
larly, the p^r electrons (yellow worldlines in fig. 

[IT] ) move approximately as in section [ll| until t + (i?—r)/c, 
i.e get the main backward boost (acceleration is maximal 
around t). Eq. (28)2 is equivalent to 


^ex ^ ^5 

or 0<{tex-l/c)vP<R 


r R; 

r R- {tex-llc)vP > 0. 


Using cylindrical coordinates {y^p^p) for X', one easily 
finds that for = Z-bA < 0 the electric force generated by 
the static charge distribution of fig. |5]b) , the associated 
potential energy mc^Ur and the left asymptotes of Ur are 

rZ2{z) 


{A,Z) = TAy^oiv) fdp 

Jo Jo 

rZ2iz 

7 dy 
Jo 


= —27re 


[p^+{y-zef] 


V{y-Ze)^+r^ 

n r^ 2 {Z) _ 

Ur{A,Z) = - J dy no(y)[V2/^+r2 


2 j 3/2 

ATTe^N{Z}. (Bl) 


(B2) 

-V(2/-^-A)2+r2] -pN{Z)A-pM{Z), 

nZ2{Z) - pZ 

Ur{-oo,Z) = - jdyn(fy)\£ZZp-y +r dyn(fy)y.{B3) 

^Jo ^ ^ Jo 

Here p^^.'K^jvfK?. Ur is continuous in (—Z,Z), since we 
have chosen U{—Z^Z) as the (A-independent) ‘additive 
constants’. Energy conservation implies 


If the left-hand side of the first line is fulfilled the surface 
electrons are expelled while the laser pulse is still entering 
the bulk and thus producing an outward force that keeps 
the LE out of Cr. Otherwise, the left-hand side of the 
second line ensures that the distance inward travelled by 


7-W,(A,Z)=7[/,Z]^,[A(/,Z),Z]=7[ee.(Z),Z]-W,(-Z,Z). 

The last equality holds only if Ze(/,Z) >0, i.e. I <fex(^)j 
the right-hand side is the electrons’ energy when expelled 
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from the bulk. This leads to the final relativistic factor 

rZ2{z) 


7/ {Z)=j{},Z} -iiN{Z)zJi,Z) + ^Jdy ndy) \^-^J[y-zJi,Z)] 


rZ2(z) 


^+r 


= ^%J,Z),Z] + ^jdyn^y-^Jy^+r^\ (B4) 

Deriving this and the identity N[Z 2 {Z)] =2N{Z) we find 
no[Z2(Z)]^ =2rLo{Z) and that, as claimed, jf{Z) is 
strictly decreasing, since d^fjdZ is negative-definite: 

% = - l^^o[Z2{Z)] 


-yndZ)zll,Z)-yN{Z) ^ f-{l,Z) 


^ dzeilZ) _ dj 
dZ 

-ym{Z) \AZ 2 iZ) - ze(i 2) ] 2+r2 - Z 2 iZ) + zjl, Z) 


For the step-shaped initial density, setting 0 = A(/,Z)—Z, 

Ur{A,Z)=f [(A-Z)v'(A-Z)2+r2-4Z(A + Z) 

+r^ sinh“''^^ — (A+Z)y(A+Z')2+r2 — 

+2Z2+2ZV4Z2+r2 +r^s iTi h-^g^] , _ 


Ur{-oo,Z)= ^ 


v4Z2+r2-Z+|^sinh-i^, (B6) 

7/(^) = + f {<^V^+^-4ZA(^,^+sinh-^f 

-[(/)+2Z]vl</>+22']2+r2-r2 sinh“'^| (B7) 




Z^— f^\4;Z2+r2 — ^ sinh ^— 


if l<iedZ). 


lil<iex{Z) then dz^ = Q = dzA at ^ = ^ex{Z), eq. ( |B5[ ) 
reduces to d'^j/dZ = M[Z — \/4Z2+r2], and (32) to 


(B5) 


l/i/(7) =MZ^[v4z2+r2-Z]^^ 


zhY 


(B8) 
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